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Combining Equations (15) and (12); 


lim N(0,t) - TGe'^ Rt-t/i 
I / L“* 


a result quoted in the Mterature (Reference 4). Evidently Equation (16) is 
valid if the diffusion length I, is small compared to the thickness of the base 
region I of the solar cell. On the oth.r hand. Equation (11), being exact within 
the limitations of the underlying theory. Is valid for all ratios '•/L. 

An alternative form to Equation (11) can be found if tlie hyperbolic functions 
of the integrand of Equation (9) are expanded in powers of exp (-2/uv /L) . The 
first term of such a series would only yield Equation (16) upon Integration, and 
subsequent terms would give corrections corresponding to situations in which 
4/1. becomes progressively smaller. The complete series must, however, become 
identical with Equation (11) upon integration. Expanding the integrand of 
Equation (9) in the manner proposed gives 


N(0, t) = T(; 


_-t/i 

'irfl 


iH-b 


du e 


ut/i 


-i'.+b 


r- r- 22 

/u(/u + iL) u - X L 


- i>’-/ut /L 


t [(■-"')■ -.ti: 

n=0 n=l 


) 


(17) 


The first term within the braces of Equation (17) gives Equation (16), as 
has been demonstrated. Since the infinite sums of Equation (17) converge 
absolutely on the integration path, summation and integration may be inter- 
changed. One is then left with a series of integrals that can be handled in the 
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ABSTRACT 


When the diffusion length of minority carriers becomes comparable with or 
larger than the thickness of a p-n junction solar cell, the characteristic decay 
of the photon-generated voltage results from a mixture of contributions with dif- 
ferent time constants. The minority carrier recombination lifetime t and the 
time constant where I is essentially the thickness of the cell and D the 

minority carrier diffusion length, determine the signal as a function of time. 

It is shown that for ordinary solar cells (n'*'-p junctions), particularly when the 
diffusion length L of the minority carriers is larger than the cell thickness 
the excess carrier density decays according to exp (-t/x-Tr^Dt/AH^) , r being the 
lifetime. Therefore, t can be readily determined by the photovolta, decay (PVD) 
method once D and f, are known. 

An analysis of this matter was published recently in the Journal of Applied 
Physics. This report offers details of its mathematical development. 
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ANALYSIS 


The author recently gave an account of the photovoltage decay (PVD) method 
for minority carrier lifetime measurements ('Reference 1) in which he promised 
to furnish the details of the mathematical analysis in a separate report. This 
report fulfills that promise. 

Since several cases (ordinary solar cells, BSF solar cells, etc.) have been 
treated In Reference 1, and since the mathematical formalism is similar for all 
cases considered, only one specific case will be considered in detail here. It 
consists of the determination of N 2 i( 0 ,t), the excess electron density at 
X = 0 (the edge of the depletion layer of the p-type base) as a function of time 
The subscript 21 identifies the mode of illumination (2 = excitation by a sliort 
flash) and the cell type (1 = an ordinary n"*'-p junction solar cell). The solar 
cell is subjected to a short pulse of duration T at t = 0; the ensuing decay of 
the excess minority carriers given by (0,t) is considered. For details, 

see Reference 1. 

Mathematically, the problem of the determination of N 2 j^ (0,t) is expressed 
thus: N 2 j (x,t) satisfies the diffusion equation in the quasi-neutral base 


^21 ■ ^21 " ’^21 ■ D~^Ge“‘'’‘T6(t) (1) 

where 

L = electron diffusion length 
D = electron diffusion constant 

-3 -1 

G = electron-hole pair generatic.i rate in cm s 
a = light absorption coefficient 

and a prime indicates differentiation with respect to x, a dot indicates 
di'^ferentlation with respect to time. 

The boundary conditions are given by 


N^^(0,t) = 0, N2j^(e,t) = 0, N2 ^(x,0) = 0 (2) 
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The first condition expresses zero current flow (open circuit); the second 
implies a perfect ohmic contact at the back surface («. being the width of the 
base) and the last, an initial condition, presupposes that the device Is 
quiescent just before it is illuminated by the pulse TG(t).* 

A Laplace transformation of Equations (1) and (2) defined by 


n(x,s) 


e ’^*'N2^(x,t)dt 


(3) 


gives the Llowing diffusion equation for the electron density in x and s 
snace: 


n" - i?‘^n 


-D 


-1 


TGc 


-nx 


(4) 


where 


- 1/2 

= L(1 + [s) (5) 

with T as the lifetime of the minority carriers. In deriving Equation (4) from 
Equation (1), use has been made of the connection between the diffusion constant 
L and the lifetime i, e.g., L = (Dt)^/^, and the rules concerning Laplace trans- 
formations (Reference 2). The solution to Equation (4) that satisfies the boun- 
dary conditions (2) may be obtained by standard methods and reads: 

n(x,s) = Ae + A[ cosh( f / 4/' ) ] {at, sinli[ (x - < )/S/’] - e cosh(x///')j (6) 


with 


A = TCid IS 


/l.h' 


l7) 


6(t) is the Dirac 6-function and T is again small comitared with the minority 
carrier recombination lifetime i or other relevant time constants. Tlie pulse 
illuminates the cell at t = i with lim / = 0 and tlie argument of the S-function 
in Equation (1) should read (t - f.) in order to make Equation (1) consistent 
with the last expression of Equation (2). 
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The electron density at the edge of the depletion layer x = 0 Is now found from 
Equation (6) by means of a back transformation (Reference 2), thus: 


N^^CO.t) 


'l«H-a 


27t1 


e^^n(0 , s)ds 


-l“4-a 


( 8 ) 


Substituting u = 1 + is and Inserting Equation (6) into this expression yields 
the following integral for N ->2 = N (omitting the subscript 21 for simplicity): 


N(0,t) 


’ i>»+b 


l(>e 7T— 
2 tt 1 - 


du — - 


Ut/l 


-l-^b 


2 2 
u a L 


e 

CO sh(v^£/L) 


tanh(»4i£/L) 

/u 


(9) 


The path of integration in the complex u-plane is taken to the right of all 
singularities of the integrand, as indicated. The integrand of Equation (9) is a 
mjromorphic function of u and has simple poles at 

u = -(n + 1/2)^ (10) 


for integer values =0, 11, 12, 13... the contour can be swung into the left- 
hand plane (Reference 2) and the residues may be picked up one by one, leading 
to the infinite series 


N(0,t) = TGe 


-t/i 


\ ^ - f (-l)'^Ti(n 

.C 4. (, 


-ut 


+ l/2)e ,2, , , ,„,2 

— ~ exp[-i! (n + 1/2) Dt/i 


n=0 


( I » ) ~ + ^ (n + 1 /2 ) ‘ 


( 11 ) 


quoted as Equation (14) of Reference 1. 


At first glance, a pole might be expected at u = n L , but its residue being 
zero, the integrand is regular at u = 
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If we let £/L go tc nfinlty in Equation (9) we obtain 


i 


lim ^(0,t) 


TCe 


-t/x 



ut/x 

e 



( 12 ) 


The value of the integral can be found in a table of integral transforms such as 
that given by Doetsch (Reference 2) or it can be evaluated directly without undue 
labor. Here tiie latter cjtion is taken; a number of similar integrals to be 
encountered below need not then be treated as extensively- 

Thfc integrand if Equation (12) has branch points at 0 and If we take the 
cut along the negative real axis and swing the contour so that the path of inte- 
gration lie-; pa/dlle"! and very close to the negative real axis in the complex 
u-plann, ir suet t iianner that the part of the path above the cut lies on the 
first Ri; mann s'leet and tlie part below the cut lies on the second Rlemann sheet 
of the double -V. lued integrand, '.e obtain 


•i<»+b 


I = 


ut/x ut/x 

^ — du = lim I du 


, , ,, i 4 j(/u + (iL) e=0 / , -/u(-/u + ah) 

/„i«>fb »/-«>- it; 


-'»+lc 


+ lim 
f+0 


ut/x 

du 


'ic 


/u(fu + aL) 


(13) 


with the substitution u = -x, this expression becomes 


I 


2iuL 



-xt/x . 

e _ 

fx(a^L^ 4- x) 



2 , 
,,“y t/ 1 


2 2 2 
y + X L 


(14) 


The last integral can be found in a table of integrals (Reference 1) with the 
result that 


I 


2TTie'^ Erfc(d/bt) 


(15) 


The complementary error function Erfc(x) is defined by Erfc(x) 



2 

- ri 

£• d n , 


4 


Combining Equations (15) and (12); 


11m N(0,t) = TGe'^ Rt-t/i ( 16 ) 

i / 


a result quoted In the Mterature (Reference 4). Evidently Equation (16) is 
valid If the diffusion ^ength L is small compared to the thickness of the base 
region I of the solar cell. On the oth-.r hand, Equation (11), being exact within 
the limitations of the underlying theory, is valid for all ratios ^ /L. 

An alternative form to Equation (11) can be found if the hyperbolic functions 
of the integrand of Equation (9) are expanded in powers of exp (-2/uv /L) . The 
first term of such a series would only yield Equation (16) upon integration, and 
subsequent terms would give corrections corresponding to situations In which 
l/L becomes progressively smaller. The complete series must, however, become 
Identical with Equation (11) upon integration. Expanding the Integrand of 
Equation (9) in the manner proposed gives 


-t/i 

N(0,t) = TG -j- 


i-bb 


du e 


ut/il 1 


-h’+b 


/u(/u + ll.) U - i L" 


i > - /uf /L 


n=() n»l 


I 

J) 


(17) 


The first term within the braces of Equation (17) gives Equation (16), as 
has been demonstrated. Since the infinite sums of Equation (17) converge 
absolutely on the integration path, summation and integration may be inter- 
changed. One is then left with a series of integrals that can be handled In the 
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saiTC manner as the integral of Equation (13). The result of the Integrations Is 
given by 


N(0,t) - TGe 


i ^ 

“ 4 rtc(«/D?) - ^ <-l)"P"“'Errc( 


«/dc + (2n + ])P/2/nt 


( 


- e'^^"'^^^“‘’^Erfc(a/Dt - (2n + l)f/2/Dt) 


+ (-l)”^e^'^”^Erfc (a*^ + n//bt) + e 

n*l 


X 


Erf c (ai'^ 


n//Dt)j| 


(IH) 


Equation (18) is the alternative expression to Equation (11) and is in fart 
Identical with it. First, since Erfc(O) = 1, Erf c (• ) “ 0 and F,rfr(-' ) = 2, the 
limit of t ■ 0 yields for Equation (18) 


11m N(0,t) 

t-0 


TG 


1 + 2e 


■2.xl 


OP 

n»0 




TG . 


(19) 


since the sums cancel. This result is Identical with the limit obtained from 
Equation (11). Furthermore, in the limit of very large absorption (it ■ 1), 
it is possible to use the asymptotic app'‘oximat Ion for Erfc 


2 

-X 

Erfc(x) = (20) 

</v-x 
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Under the aHsumptlon that a/Ot '^nS/Zot, it can he ahown that Equation (18) 
becomes 


Urn N(0,t) 


-t/T 

/■nDt 


( f m ^ 

1 2 ^ (- 1 )" 

< n-1 ! 


( 21 ) 


But this expression is related to a theta function (Reference 5), 


\(l/2, i Dt/.^) = — '^^-(l + 2y^ (-1)" e"" ' 


( 22 ) 


so th.it Kijuation (21) becomes 


1 im N(0,t) 

I* * 


[Gv 


-t/ 


ai 1 


(1/2, lrl)t/t ) 


(23) 


Kurthi-rmore (K«-terence 3;, 


(v,:) = 2 ^ ' (-1)" sin[(2n H l)iiy 

n-O 


(24) 


that finallv 


n*() 


(25) 


an expression that is identical with Equation (’I), provided that at ■■ 1. Since 
both series. Equation (11) and Equation (21), identical, the que.stio'’ ari.ses: 

which is to he preferred? As an example, take a “ 10^/cm, t “ 10"^/cm and 

I) = 3b t m-/sec, correspond in^t to values for Si at a duping level of 

^ l()l'’/cm^ and a wavelength of 0.55 um (maximum solar intensity). If t is 
measuri'd In ..sec, the arguments of the complementary error functions of 
Equation (18) become bOrt ’3/2/t and f or t 1 Msec, the approximation (25) 

already is verv good. li i is small, however, the photon energy being near the 

band gap, manv terms of Equation (18) ,.iust be used for longer times; for the 
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series (11), only the first term need be retained. Therefore, Equation (11) la 
to be pieferred. It was used in Reference 1. 

Calculations similar to those made on these pages for the other cases dis- 
cussed in Reference 1 lead to analogous results. 
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